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Because of recent interest in bats as reservoirs of 
emerging diseases, we investigated the presence of viruses 
in Mystacina tuberculata bats in New Zealand. A novel 
alphacoronavirus sequence was detected in guano from 
roosts of M. tuberculata bats in pristine indigenous forest on 
a remote offshore island (Codfish Island). 


He settlement in New Zealand is relatively recent 
compared with that of many countries and extends back 
~800 years to when early Polynesian explorers first arrived 
(1). A wide variety of exotic flora and fauna have since been 
introduced, and major changes to the landscape and ecology 
have occurred, particularly once European settlers arrived 
150 years ago. However, large areas of New Zealand that are 
representative of a prehuman state still remain, in particu- 
lar offshore, islands such as Whenua hou (Codfish Island), 
which is situated off the southern coast of New Zealand. Be- 
fore human settlement, only 3 species of terrestrial mammals 
were present, all of which were bats. Why nonvolant mam- 
mals have been absent is unknown, but this absence was 
probably caused by major extinction events and subsequent 
geographic isolation that prevented recolonization (2,3). 
Only 2 bat species remain in New Zealand and both 
are considered vulnerable (www.iucnredlist.org/). The 
long-tailed bat (Chalinolobus tuberculatus) belongs to the 
family Vespertilionidae and is believed to have arrived 
from Australia ~1 million years ago (4). The lesser short- 
tailed bat (Mystacina tuberculata) is the sole surviving 
member of the family Mystacinidae in New Zealand (5). M. 
tuberculata bats are believed to have diverged from other 
bat species ~20 million years ago (6) and have lived in 
geographic isolation from other bat species until the arrival 
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of C. tuberculatus bats (3). M. tuberculata bats also exhibit 
one of the widest feeding ranges of any bat species and 
occupy a niche similar to rats and mice because these bats 
can walk on the ground by using their wings (7). 

Little is known about the microorganisms present in 
these bats. One study found no evidence of pathogenic 
bacteria or lyssaviruses in them, but reported a Sarcocystis 
sp. (8). Given the recent intense interest in bats as a reservoir 
of emerging diseases and the advent of high-throughput 
sequencing as a virus discovery tool, we investigated the 
presence of viruses in M. tuberculata bats. 


The Study 

Four bat guano samples were collected from known 
roost sites on the remote offshore island of Whenua hou 
(Codfish Island) (46°47'S, 167°38’E), which is situated at 
the southern coast of New Zealand. This small island is 
heavily forested and largely unmodified by humans. It has 
special conservation status for the protection of endangered 
species, and public access is not permitted (9). 

Bat guano was held at 4°C during transport (<48 
h), resuspended in 2 mL phosphate-buffered saline, and 
subjected to centrifugation at 6,000 x g for 5 min. RNA 
was extracted from 400 uL of supernatant by using the 
iPrep PureLink Virus Kit (Life Technologies, Carlsbad, 
CA, USA) and eluted into 50 uL reverse transcription 
PCR molecular-grade water (Ambion, Austin, TX, USA). 
Metagenomic sequencing was then conducted for 1 of the 
samples by using an Illumina MiSeq Instrument (New 
Zealand Genomics Ltd., Massey Genome Service, Massey 
University, Palmerston North, New Zealand) after a series 
of steps involving DNase I treatment, reverse transcription, 
multiple displacement amplification (QIAGEN, Valencia, 
CA, USA), and Illumina TruSeq library preparation (New 
Zealand Genomics, Ltd.). 

A total of 10,749,878 paired-end sequence reads 
of 250 bp were generated and assembled into contigs by 
using Velvet 1.2.07 (10). Assembled contigs were searched 
for viral sequence by comparison to the nonredundant 
nucleotide database in Genbank (downloaded March 2013) 
by using the nucleotide basic local alignment search tool 
(National Center for Biotechnology Information, Bethesda, 
MD, USA). Forty-six contigs showed similarity to known 
genus Alphacoronavirus sequences (Table). (Raw sequence 
data are available on request from the authors.) 

Virus sequence from the genus A/phacoronavirus was 
confirmed in all 4 original guano samples by using a specific 
reverse transcription PCR based on the metagenomic 
data specific for 582 bp of the RNA-dependent RNA 
polymerase (RdRp) gene (GenBank accession nos. 
KF515987-KF515990). The Rdrp sequence was identical 
in all 4 guano samples, and the closest relative was bat 
coronavirus HKU8 (GenBank accession no. DQ249228; 


697 


DISPATCHES 


Table. Summary of BLASTn output for contigs from bat guano that show identity to alphacoronaviruses, New Zealand* 


Highest scoring BLASTn hit recorded against 


No. contigs matching Range of contig lengths, alphacoronaviruses 
Alphacoronavirus gene alphacoronaviruses bp e-value Nucleotide identity (%)t 
Open reading frame 1ab 33 182-1054 8x10 581/828 (77) 
Spike protein 5 580-1629 1x10 362/551 (66) 
Matrix 4 251-840 2x10 532/746 (71) 
Nucleocapsid 4 536-890 6x 10° 79/109 (72) 


*BLASTn, nucleotide basic local alignment search tool 


(http://blast.ncbi.nlm.nih.gov/Blast.cgi? PROGRAM=blastn&PAGE_TYPE=BlastSearch&LINK_LOC=blasthome). 


tNo. nucleotides identical between query and database sequence. 


79% nt identity, 426/542). Phylogenetic analysis was 
performed for Rdrp (Figure 1) and for the spike protein, as 
derived from metagenomic data (GenBank accession no. 
KF575176) (Figure 2). We propose that these data support 
identification of a new alphacoronavirus, which has been 
designated as Mystacina bat CoV/New Zealand/2013. 


AY994055 Feline infectious peritonitis virus 
JN856008 Canine coronavirus strain A76 
JQ408981 Feline infectious peritonitis virus strain DF-2 
NC 002306 Feline infectious peritonitis virus 
JQ408980 Feline coronavirus strain DF-2 R3i 
94 AJ271965 Transmissible gastroenteritis virus 
| DQ811786 TGEV Miller M60 
9766 0Q811785 TGEV Miller M6 
DQ811787 PRCV ISU-1 
98 KC175340 Canine coronavirus strain K378 
FJ938051 Feline coronavirus RM 
HM245925 Mink coronavirus strain WD1127 
93 HM245926 Mink coronavirus strain WD1133 
NC 009988-Rh Bat coronavirus HKU2 
52 AF304460 Human coronavirus 229E 
99° JQ410000 Alpaca respiratory coronavirus isolate CA08-1/2008 
NC 010438 Bat coronavirus HKU8 
gg | ©U420137 Bat coronavirus 1B strain AFCD307 
85 EU420138 Bat coronavirus 1A strain AFCD62 
69 JQ989271 Rousettus bat coronavirus HKU10 isolate 175A 


109 JQ989272 Hipposideros bat coronavirus HKU10 isolate TLC1343A 
JQ989268 Hipposideros bat coronavirus HKU 10 isolate TLC1310A 


JQ765566 Human coronavirus NL63 strain NL63/DEN/2008/16 
JQ765568 Human coronavirus NL63 strain NL63/DEN/2005/193 
JX524171 Human coronavirus NL63 isolate CBJ123 
JQ765564 Human coronavirus NL63 strain NL63/DEN/2009/14 
JQ765570 Human coronavirus NL63 strain NL63/DEN/2005/235 
NC 009657 Scotophilus bat coronavirus 512 
109. 1 NC 003436 Porcine epidemic diarrhea virus 
9 JQ023162 Porcine epidemic diarrhea virus strain attenuated DR13 
JX088695 Porcine epidemic diarrhea virus strain GD-B 
JX560761 Porcine epidemic diarrhea virus strain SD-M 
KF515987 Mystacina bat CoV/New Zealand/1/2013 
KF515988 Mystacina bat CoV/New Zealand/3/2013 
87 KF515989 Mystacina bat CoVNew Zealand/2/2013 
KF515990 Mystacina bat CoV/New Zealand/4/2013 
EF065513-Ro-batCov Bat coronavirus HKU9-1 
196 AY772062 SARS coronavirus WH20 
59 DQ22305-Bat SARS coronavirus HKU3-1 
EF065509-PiBatCoV Bat coronavirus HKUS-1 
DQ415908-rdrp Human coronavirus HKU1 strain N11 genotype A 
gg | AY585228-HCov Human coronavirus OC43 strain ATCC VR-759 
57, AF220295-BCoV Bovine coronavirus strain Quebec 
53 EF 424623 Giraffe coronavirus US/OH3/2003 
DQ646405-IVB Infectious bronchitis virus isolate TW2575/98 
95 EU111742-SW1 Coronavirus SW1 
54 FJ376619-BuCoV Bulbul coronavirus HKU11-934 
FJ376621 Thrush coronavirus HKU12-600 
100 FJ376622-MunCoV Munia coronavirus HKU13-3514 
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85 JQ065042-PorCov Porcine coronavirus HKU15 strain HKU15-44 4 


0.05 


Conclusions 

The discovery of unknown coronaviruses provides 
information for a model of coronavirus evolution (//) 
and contributes to understanding the process of disease 
emergence, as in detection of Middle East respiratory 
syndrome coronavirus (/2). The  alphacoronavirus 


7 Figure 1. Phylogenetic tree showing 
genetic relatedness of RNA-dependent 
RNA polymerase amino acid sequences 
for Mystacina sp. bat coronavirus (CoV)/ 
New Zealand/2013 (shown in boldface) 
with those of known coronaviruses. 
Evolutionary history was inferred for 
183 informative amino acid sites by 
using the maximum-likelinood method 
based on the Whilan and Goldman 
model with gamma distribution in MEGA 
5.05 software (www.megasoftware. 
net). Bootstrap values are calculated 
from 1,000 trees (only bootstrap values 
>50% are shown). Scale bar indicates 
- a-CoV nucleotide substitutions per site. TGEV, 
transmissible gastroenteritis CoV; PRCV, 
porcine respiratory CoV; SARS, severe 
acute respiratory syndrome. 
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NC 009988-Rh Bat coronavirus HKU2 + a-CoV Figure 2. Phylogenetic tree 
DQ646405-IBV Infectious bronchitis virus isolate TW2575/98 “L yCov showing genetic relatedness 
EU111742-SW1 Coronavirus SW1 J of spike protein amino acid 
92 DQ811785 TGEV Miller M6 7] sequence for Mystacina sp. 
re DQ811786 TGEV Miller M60 bat coronavirus (CoV)/New 
pened lie tae acl ae Zealand/2013_ (shown in 
AJ271965 ‘Trahemiaaiie gastroenteritis ius boldface) with those of known 
JQ408980 Feline coronavirus strain DF-2 R3i F 
ae JQ408981 Feline infectious peritonitis virus strain DF-2 Se eae Evolutionary 
99 AY994055 Feline infectious peritonitis virus history was inferred for 492 
65 NC 002306 Feline infectious peritonitis virus informative amino acid sites by 
KC175340 Canine coronavirus strain K378 using the maximum-likelihood 
99) | HM245925 Mink coronavirus strain WD1127 method based on the Whilan 
95 99 HM245926 Mink coronavirus strain WD1133 and Goldman + F_ model 
JN856008 Canine coronavirus strain A76 with gamma distribution and 
FJ938051 Feline coronavirus ep invariant sites in MEGA 5.05 
NC 009657 Scotophilus bat coronavirus 512 software (www.megasoftware. 
88 JQ023162 Porcine epidemic diarrhea virus strain attenuated DR13 
98 _ JX560761 Porcine epidemic diarthea virus strain SD-M net). Bootstrap values are 
99; NC 003436 Porcine epidemic diarrhea virus - a-CoV calculated from 1,000 trees 
65 JX088695 Porcine epidemic diarrhea virus strain GD-B (only bootstrap values >50% 
7 55 NC 010438 Bat coronavirus HKU8 are shown). Scale bar indicates 
KF575176 Mystacina bat CoV/New Zealand/3/2013 nucleotide substitutions per 
86 100 EU420137 Bat coronavirus 1B strain AFCD307 site. TGEV, transmissible 
Set _ ee 1" agp eee gastroenteritis CoV; PRCV, 
‘ousettus bat coronavirus isolate * é 
94 JQ989268 Hipposideros bat coronavirus HKU10 isolate TLC1310A porcine —_ respiratory CoV; 
100 JQ989272 Hipposideros bat coronavirus HKU10 isolate TLC1343A SARS, severe acute respiratory 
50 98 AF304460 Human coronavirus 229E syndrome. 
JQ410000 Alpaca respiratory coronavirus isolate CA08-1/2008 
73 JQ765564 Human coronavirus NL63 strain NL63/DEN/2009/14 
a JQ765570 Human coronavirus NL63 strain NL63/DEN/2005/235 
99 JQ765568 Human coronavirus NL63 strain NL63/DEN/2005/193 
80 JQ765566 Human coronavirus NL63 strain NL63/DEN/2008/16 
88 JX524171 Human coronavirus NL63 isolate CBJ123 J 
FJ376621-ThCoV Thrush coronavirus HKU12-600 > 
FJ376619-BuCoV Bulbul coronavirus HKU11-934 L §-CoV 
99 FJ376622-MunCov Munia coronavirus HKU13-3514 
66 JQ065042-PorCoV Porcine coronavirus HKU15 strain HKU15-44 J 
EF065513-Ro-batCoV Bat coronavirus HKU9-1 “| 
99 EF065509-PiBatCoV Bat coronavirus HKU5-1 
= 99° AY772062 SARS coronavirus WH20 
DQ022305-Rh Bat SARS coronavirus HKU3-1 & B-C oV 
88 DQ415908-HCoV Human coronavirus HKU1 strain N11 genotype A 
100 AY585228-HCoV Human coronavirus OC43 strain ATCC VR-759 
96|| AF220295-BCoV Bovine coronavirus strain Quebec 
58 EF424623 Giraffe coronavirus US/OH3/2003 =) 


0.5 


identified in this study from M. tuberculata bat guano 


is unique in respect to the extreme geographic and 
evolutionary isolation of the host bat species, which along 
with C. tuberculatus bats, has been separated from all other 
mammalian species for ~1 million years. 

The current estimate for a common ancestor for all 
coronaviruses is 8,100 Bc (/3). To be consistent with this 
estimate, Mystacina bat coronavirus would need to have 
been introduced to bats on Whenua hou within the past 800 
years since humans first arrived on this island (given that 
the island had no other mammals before this time) (//) or is 
extant to modern alphacoronavirus phylogenetic radiation 
(genesis and expansion). Apart from humans, only 2 
other terrestrial mammals have ever inhabited Whenua 
hou: the brushtail possum (Trichosurus vulpecula) and 
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the Polynesian rat (Rattus exulans), both of which were 
eliminated from the island in the late 1980s (//); neither 
mammal has been reported as a host of alphacoronaviruses. 
Members of the genus Alphacoronavirus infect only 
mammals. Thus, an avian origin for this virus is unlikely. 
An alternative theory of an ancient origin for all 
coronaviruses has recently been proposed that involves 
an alternative evolutionary molecular clock analysis, 
which places the most recent common ancestor many 
millions of years ago (/4). The discovery of Mystacina 
bat CoV/New Zealand/2013 virus could lend support 
to such a theory; despite potentially millions of years of 
isolation, it has diverged relatively little from other extant 
alphacoronaviruses, as shown by the close relationship of 
the Rdrp and spike protein genes to those of other extant 
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alphacoronaviruses (Figures 1, 2). An expanded survey for 
Mystacina bat coronavirus in mammals in New Zealand 
and subsequent characterization of viral genomes would 
provide further insights into the origin of coronaviruses. 
No instances of human zoonotic disease from New 
Zealand bat species have been reported. Only a small number 
of conservation staff handle these bats, these staff use 
standard personal protective equipment and work practice. 
Staff are also offered prophylactic rabies vaccination as a 
precautionary measure, even though New Zealand is free 
from rabies. The genus Alphacoronavirus includes several 
human and animal pathogens, but on the basis of phylogenetic 
data in this study, it is not possible to estimate the risk posed 
by Mystacina bat coronavirus to human or animal health. 
New Zealand is not considered a hot spot for emerging 
infectious diseases. This country is free from many human 
and animal diseases, such as rabies and foot-and-mouth 
disease, and infections with human arboviruses because of 
recent colonization by humans and strict biosecurity border 
controls (/5). Thus, indigenous wildlife in New Zealand 
has generally been viewed as an almost sterile and unique 
biosphere. Given detection of this coronavirus, more 
thorough characterization of the ecology of viruses and other 
microorganisms in native wildlife should be considered to 
fulfill conservation needs and further safeguard human and 
domestic animal health against cross-species transmission. 
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